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ABSTRACT 

We argue that the analog nature of quantum computing makes the usual design approach of constructing complicated 
logical operations from many simple gates inappropriate. Instead, we propose to design multi-spin quantum gates 
in which the input and output two-state systems (spins) are not necessarily identical. We outline the design criteria 
for such devices and then review recent results for single-unit Hamiltonians that accomplish the NOT and XOR 
functions. 
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1. INTRODUCTION 

One of the great challenges of the physics of nanoscale systems has been the design of atomic-size devices oper- 
ating in a quantum-coherent fashion. Dimensions of semiconductor computer components will soon reachOI about 
0.25 /im = 2500 A, which is well above the atomic sizes at which quantum-mechanical effects are important. However, 
it is generally expected that as the miniaturization continues, atomic dimensions will be reached. This article con- 
cerns with quantum computin^Jj2_ nanoscale devices that perform logical operations while maintaining quantum 
coherence. Some early studiesaau'ElEI considered how quantum mechanics affects the foundations of computer science; 
issues such as^lipaiuyinnis rmiibjissical cninpulatinn dui2±D qiia|nt]iin. fl.uctu;^ etc., have been raised. A more recent 
developmentymlaMme y ymBB'B has been to utilize the quantum nature 

of components of atomic dimensions for more efficient computations involving quantum-coherent evolution. 

Quantum computing has attracted a lot of interest recently owing to several 
faster t,>i{jn pl^jsi^jr?)! pnpnnnting- new fast quantum algorithms have been proposed 
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have been explored. Ideas on how to combine the 
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and some Hamiltonians for gate operation 
simplest quantum gates have been put forth 
simplest quantum-gate functions have been recently realizec 

There remain, however, many conceptual difficulties with quantum computing.l 
quantum evolution implies that the time scale At of the operation of quantum logic gates must be built into the 
Hamiltonian. As a result, all the proposals available to date assume that computation will be externally timed, i.e., 
interactions will be switched on and ojf, for instance, by laser radiation. 

This means that if logical operations are constructed from one or few simple universal gates, then each such gate 
will have to be precisely controlled from outside. In ordinary classical (i.e., macroscopic, irreversible) computing, the 
NAND gate is an example of a universal gate. From it complicated logical operations can be constructed. In the 
classical case, however, it is the internal relaxation processes in the basic gate(s) that determine the time scale of their 
operation (equilibration) At. We consider it extremely unlikely that one would ever be able to control externally, in 
a coordinated fashion, millions of simple reversible quantum gates in order to operate a macroscopic computer. 

Furthermore, quantum computers are naturally ana^ofyii in their operation. Indeed, in order to use the power 
of quantum interference (superposition of states), one has to allow any linear combination of the basis qubit states 
|1) and |0). Analog errors are difficult to correct. By analog errors we mean those minor variations in the input 
and output variables which cannot on their own be identified as erroneous in an analog device because its operation 
involves continuous values of variables (so that the fluctuated values are as legal as the original ones). By noise errors 
we term those that result from single-event problems with device operation, or from external influences (including 



decoherence in the quantum case), or from other failures in operation. All errors in a digital device (i.e., deviations 
from discrete values) can be systematically decreased or eliminated in each step of a calculation. Similarly, the noise 
errors in analog devices can be corrected or decreased. 

However, the analog errors cannot be corrected. Consider a state a|l) + f3\Q) and a nearby state + /3'|0) + 
Cj where a' is close to a, (3' is close to /3, while Q are small. The latter terms represent admixture of quantum 
states |j) other than the two qubit states. Both states are equally legal as input and output quantum states. We could 
restrict input or output to a vicinity of certain states, for instance, the basis states |1) and |0), thus moving towards 
digitalization. However, we then loose the quantum-interference property. Another important effect: decoherence, 
that would require a density matrix description, falls in the noise-error category. 

Modern error-correction technique can handle the noise errors but not the analog errors. To 

illustrate, consider this quoteea from the article entitled Quantum Error Correction for Communication: "To achieve 
this the sender can add two qubits, initially both in state |0), to the original qubit and then perform an encoding 
unitary transformation. . ." . The problem here is that the states actually encountered in the system during error 
correction are not available as basis qubit states (such as |0)) with infinite precision. Typically, by qubits we mean 
a set of two orthogonal quantum states selected from the energy eigenstates of the system. Even assuming that the 
thermal noise can be reduced at low temperatures to make the ground state sufficiently long-lived, the excited states 
of any system, especially if it is a part of a macroscopic computer, will not be defined sharply enough to provide 
ideal stationary states |1) and |0). External interactions, spontaneous emissions, etc., will generate both noise- and 
analog-errors in the basis states, i.e., the actual state (disregarding decoherence) will be a\l) + 13\0) Cj|j)i with 

a ~ 1 and /3, Q — 0, instead of the ideal |1) which is an eigenstate of an ideal, isolated-system Hamiltonian. 

Furthermore, analog errors will be magnified when separate simple-gate operations are combined to yield a 
complex logical function. Thus, the conventional picture of a quantum computer is unrealistic: it assumes a multitude 
of simple-gate units each being externally controlled by laser beams (one needs a lot of graduate students for that!). 
Such computers will magnify analog errors which cannot be corrected in principle because the error state is as legal 
as the original state. 

In this work we therefore adopt a view typical of the analog-computer approach, of designing the computer as 
a single unit performing in one shot a complex logical task instead of a chain of simple gate tasks. This approach 
will not repair all the ailments outlined earlier. For instance, the computer as a whole will still be subject to analog 
errors. However, these will not be magnified by proliferation of sub-steps each of which must be exactly controlled. 

In fact, we consider it likely that technological advances might first allow design and manufacturing of limited-size 
units, based on several tens of atomic two-level systems, operating in a quantum-coherent fashion over a sufficiently 
large time interval to function as parts of a larger classical (dissipative) computer which will not maintain a quantum- 
coherent operation over its macroscopic dimensions. We would like these to function as single analog units rather 
then being composed of many gates. 

The outline of this review is as follows. In Section 2 we continue our discussion of the design of quantum gates. 
In Section 3 we review known results for the simplest NOT gate mainly to set the notation and nomenclature. 
A more complicated, two-spin NOT gate is studied in Section 4. Section 5 addresses the time-dependence of the 
Hamiltonians. Finally, Sections 6, 7, 8 review results for a three-spin XOR gate. 



2. DESIGN CONSIDERATIONS FOR MULTI-SPIN QUANTUM GATES 

In order to make connections with the classical computer-circuitry design and identify, at least initially, which multi- 
qubit systems are of interest, we propose to consider spatially extended multi-spin quantum gates with input and 
output qubits possibly different. The reason for emphasizing this property is that multi-spin devices will have spatial 
extent. The interactions that feed the input need not be identical to those interactions/measurements that read off 
the output. Furthermore, for systems with short-range interactions one can only access the boundary spins in a large 
cluster. Thus we may use only part of the spins to specify the input and another subset to contain the output. The 
two sets may be identical, partially overlapping, or nonover lapping. 

Reversibility of coherent quantum evolution makes the distinction between the input and output less important 
than in irreversible computer components. However, we consider the notion of separate (or at least not necessarily 



fully identical) input and output useful within our general goals: to learn what kind of interactions are involved and 
to consider also units that might be connected to/as in classical (dissipative) computer devices. 

Our goal is to be able to design interaction parameters, presumably by numerical simulations, to have such gates 
perform useful Boolean operations. This is not an easy task. Actually, it must be broken into several steps. First, 
we must identify those interactions which can be realized in solid state or other experimental arrangements. As 
examples below illustrate for several simplest gates the form of the interaction Hamiltonians is quite unusual by the 
solid-state standards. 

Secondly, we expect interactions to be short-range and two-particle (two-spin) when several two-state systems 
(termed qubits, spins) are involved. 

Thirdly, incorporating designed coherent computational units in a larger classical computer will require a whole 
new branch of computer engineering because the built-in Boolean functions will be complicated as compared to the 
conventional NOT, AND, OR, NAND, etc., to which computer designers are accustomed. Furthermore, the rules 
of their interconnection with each other and with the rest of the classical computer will be different from today's 
devices. 

Our initial studies have been analytical. In the future we foresee numerical studies of systems of order 20 to 
25 twq-stHt|P .(spin) atomic components with variable general-parameter interactions. In this review we summarize 
resultaHS'LilES for interaction Hamiltonians required for operation of the NOT (Sections Spi^Land XOR (Sections 6, 
7, 8) kigic-gates. Other results presently available include Hamiltonians for certain NOTIiJnJ and controUcd-NOT 
gates,E3E3u3 and for snrnp copying processes as well as general analyses of possibility of construction of quantum 
computing systems. eIl3 

Quantum logics and the dynamics of quantum gates should be fully reversible. Implications of this property have 
biased recent literature on the quantum logic gates. Firstly, the distinction between the input and the output parts of 
the system has been blurred. A typical configuration involves a quantum-mechanical system that is "programmed" 
with the input and then after the time interval At it will be in the output state. We note that the time interval 
Ai is fully determined by the parameters of the Hamiltonian; in order to effect the quantum gate operation, the 
interaction energies associated with both the internal and external-field parts of the Hamiltonian must be of order 

n/At. 

Consideration of multi-spin quantum gates requires a large number of basis states. However, it is also useful to 
study few-spin exactly solvable systems. These provide explicit examples of what the actual interaction Hamiltonians 
should look like. A notable exactly solvable jSystem, known before the quantum-computing field became active, is 
the NOT gate operation in a two-state qubito obtained by applying a constant external magnetic field to a single 
spin. Then another field is applied, oscillating in time, in a direction perpendicular to the constant field. This 
paramagnetic resonance problem is a textbook example of time-dependent quantum-mechanical evolution. 

An accepted approach has been to consider interactions switched on only for the duration of the gate operation 
At. If the gate is actually the whole computer then one can regard the interaction as time- independent. However, for 
specific tasks in components with a limited number of basis states, it may be appropriate to view the interaction as 
controlled externally to be switched on and off. While general ideas of externally timed computation are not new,cl 
actual realizations in quantum computation with many sub-unit gates will encounter difficulties^oiitlined earlier. 
General developments for the latter type of interaction (time-independent or on/off) have includeaj'ES identification 
of unitary operators that correspond to quantum computer operation and establishment of the existence of the 
appropriate interaction Hamiltonians. 

A quantum gate performs an operation whereby the input state determines the output state after a time interval 
At. The interactions must be controlled, i.e., switched on and off, in order to have the gate operation during the 
interval At independent of the interactions with the computer parts external to the gate. This control of interaction, 
i.e., external timing of the computer operation already mentioned earlier, can be possibly accomplished by the 
external interactions while the internal interactions be reserved for the gate operation. However, we would like to 
consider multi-spin gates in order to avoid too many such controlling external influences. 

With regards to the requirement to control the interactions externally, with the time dependence given by the 
on/off protocol, we will show in Section 5 how to extend this approach to certain time-dependent interactions 
(protocols) which are more smooth than the on/off shape. 



3. THE SIMPLEST NOT GATE 



In this spf-|.iQTi^wj^ rnnsirlpr the NOT gataHs based on a single qubit. This gate has been extensively studied in the 
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literature 
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Our presentation here is only intended to set up the notation and illustrate methods useful in 



more complicated situations addressed in Sections 4, 6, 7, 8. We label by J and j the two basis states. The 
NOT gate corresponds to those interactions which, over the time interval At, accomplish the following: 



(1) 



JfJ 



Here a and /3 are arbitrary. The unitary matrix U that corresponds to this evolution is 



U = 



e 

s*" 



(2) 



(3) 



with the eigenvalues 



= e^(a+/3)/2 and U2 = -e'("+«/2 , 
while the normalized eigenvectors yield the transformation matrix T which can be used to diagonalize U: 



(4) 



T = 



1 



,i/3/2 pil3/2 



V2 



./2 



We have 



ia/2 



(5) 



Here the dagger superscript denotes Hermitian conjugation. 
We next use the relation 



(6) 



U = e-'"^'/^ (7) 
for a time-independent Hamiltonian. In the diagonal representation, it yields the energy levels: 



where A^i and N2 are arbitrary integers. The Hamiltonian is then obtained from the relation 



(8) 



(9) 



as a certain 2x2 matrix. The latter is conveniently represented is terms of the unit matrix X and the conventional 
Pauli matrices ax, Oy, a^'- 



H 



-2A^(" + ^)+A^(^^ 



N2 + - 



■kTi , 



N2-- 



a — P 



(10) 



To effect the gate operation, the interaction must be switched on for the time interval At. The constant part of 
the Hamiltonian only affects the average phase of the transformation (l)-(2). Thus this term can be disregarded. 

The nontrivial part of (10) depends on the integer N = Ni— N2 which is arbitrary, and on one arbitrary variable 
Thus we can use the Hamiltonian 

^ " S " + ^^'"^ ""^^ • ^^^^ 

For a spin-i two-state system such an interaction can be obtained by applying a magnetic field oriented in the 
xy-planc at an angle 7 with the x-axis. The strength of the field is inversely proportional to the desired time interval 
Ai, and various allowed field values arc determined by the choice of N . 

We note that during application of the external field the up and down quantum states in (l)-(2) arc not the 
eigenstates of the Hamiltonian. If the time interval At is not short enough, the energy- level splitting \Ei — -E2I cx; 
|7V — ^1 can result in spontaneous emission which is one of the sources of errors in computer operation. Generally, 
when implemented in condensed matter, the two states of the qubit may be part of a spectrum of many energy levels. 
In order to minimize the number of spontaneous transition modes, the best choice of the interaction strength would 
correspond to minimizing \Ei — i?2|, i-e., to |iV — ^ \ = 5. 

4. THE SPATIALLY EXTENDED TWO-SPIN NOT GATE 

In this section we consider a more complicated situation. Two two-state systems, input (/) and output (O), are 
involved. We will use the following self-explanatory notation for the state vector: 



(11) 



ail TT) + a2| Ti) + a3| iT) + a4| ii) 




In what follows we will omit the direct-product symbols ® when multiplying expressions with subscripts / and O. 

We seek transformations such that irrespective of the initial state of O, and provided / is initially in the up or 
down state, the final state has O in the down or up state, respectively (while the final state of / is not restricted). 
Thus, we place the logical NOT of / in O provided the initial state of / was one of the basis states corresponding 
to the classical bit values 1 and 0. The desired transformation maps any state with 03 = 04 = into a state with 
components 1 and 3 equal zero, i.e., input up yields output down. Similarly, any state with ai = a2 = should 
evolve into a state with components 2 and 4 equal zero, corresponding to input down giving output up. The general 
evolution operator must therefore be of the form 

/ [/i3 Uii 

jj ^ U21 U22 
loo C/33 Usi 
\U4i U42 

which depends on 16 real parameters. However, one can show that the requirement of unitarity, WU = 1, imposes 
8 conditions so that the number of real parameters is reduced to 8. The following parametrization covers all such 
unitary matrices: 





u 



f Q e^^smn 
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e*" cos T 



e^^ sin T 







e^^ cos ^ 


.g'C/J+'S-x) sin n 




(15) 



Here all the variables are unrestricted; we could limit Vl and T to the range [O, -l] without loss of generality. 

In order to make the calculation analytically tractable, we will restrict the number of free parameters to four by 

considering the matrix 



U 












ei/3 





g*p 































(16) 



This form has been favored for a possible analytical calculation for the following reasons. Firstly, the structure of a 
single phase-factor in each column is similar to that of the two-dimensional (single-spin) matrix encountered earlier. 
Secondly, the form (16) contains Hermitian-Z7 cases (/3 = —a, p = or tt, ^ = or tt). Therefore, the eigenvalues, 
which arc generally on the unit circle for any unitary matrix, may be positioned symmetrically with respect to the 
real axis, as functions of the parameters. Indeed, the eigenvalues of U turn out to be quite simple: 



(17) 



The (unitary) diagonalizing matrix T is 



1 

71 



^g^/3/2 

V2 

V2 I ■ 

\(,ia/2 _gia/2 g 



(18) 



The next step is to identify the energy levels. We chose the notation such that the energies Ei^2 are identical to 
(8). The other two energies are given by 



n 2nn _ 



E4 = --^s 

At 



2Trh . 
'At 



-Na 



(19) 



The Hamiltonian is then obtained as in the single-spin NOT case. It is convenient to avoid cumbersome expressions 
by expressing it in terms of the energies; the latter will be replaced by explicit expressions (8), (19) when needed. 
The resulting 4x4 matrix has been expressed in terms of the direct products involving the unit matrices and the 
Pauli matrices of the input and output two-state systems. We only report the result: 



H=-^ {2Ei + 2E2 +Ei + Ei) + -^ {E3 - E4) {a,i - a,o) + \ {^E^ + 2^2 -E^- E4) a,ia,o 

+ i (^Ei - E2) ^COS {(TxK^xO - (^yK^yo) + ^ (-E'l - -^2) ^sin " ^ {(TxK^yO + (^yK^xo) ■ (20) 

The constant part of the Hamiltonian can be changed independently of the other coupling constants and it can be 

discarded. We can also generally vary the integers A'^1,2,3.4 and the variables a, [3, p, S. The constant part is in 
fact proportional to 1/ (8) Jq- We avoid this notation and present the terms in the Hamiltonian in a more physically 
transparent form. 

The Hamiltonian in (20) has also terms linear in the Pauli matrices (in the spin components for spin systems). 
These correspond to interactions with externally applied fields which in fact must be of opposite direction for the / 



and O spins. As explained in the introduction, we try to avoid such interactions: hopefully, external fields will only 

be used for clocking of the computation, i.e., for controlling the internal interactions via some intermediary part of 
the system connecting the / and O two-state systems. Thus, we will assume that E3 = E4 so that there are no terms 
linear is the spin components. 

Among the remaining interaction terms, the term involving the 2-components in the product form a^jazO (= 
(^zi ® <^zo), has an arbitrary coefficient to be denoted —£. The terms of order two in the x and y components have 
free parameters similar to those in (11)- (12). The final expression is 



(COS7) (CTa,/(TccO - (^yl^vO) + (sin7) {axiCTyO + (^yK^xo) ■ (21) 



Here N = Ni — N2 must be an integer. In order to minimize the spread of the energies Ei and E2 we could choose 
I A'' — i| = i. Recall that we already have Es = E4. Thus the energy levels of the Hamiltonian in (21) are 

E^ = -£+?:(^-l)^ ^2 = -^-SfA^4V E,,4 = S. (22) 



At\ 2/ ' At\ 2, 

Thus degeneracy of three levels (but not all four) can be achieved by varying the parameters. 

The form of the interactions (21) is quite unusual as compared to the traditional spin-spin interactions in con- 
densed matter models. The latter usually are based on the uniaxial (Ising) interaction proportional to az<Tz, or the 
planar Xy-modcl interaction proportional to axCx + <^y<^y, or the isotropic (scalar-product) Heisenbcrg interaction. 
The spin components here are those of two different spins (not marked). The interaction (21) involves an unusually 
high degree of anisotropy in the system. The x and y components are coupled in a tensor form which presumably 
will have to be realized in a medium with well-defined directionality, possibly, a crystal. 



5. COMMENT ON TIME-DEPENDENCE OF INTERACTIONS 

The Hamiltonians considered thus far were all constant for the duration of the gate operation. We note that the 
external control of the interaction need not be limited to the time-dependence which is an abrupt on/off switching. 
Indeed, we can modify the time dependence according to 



H{t) = f{t)H, (23) 

where we use the same symbol H for both the original time- independent interaction Hamiltonian such as (21) and 
the new, time-dependent one, H{t). The latter involves the protocol function f{t). The shape of this function, 
nonzero during the operation of the gate from time t to time t + At, can be smooth. 

For Hamiltonians involving externally applied fields, such as (12), it may be important to have a constant plus an 
oscillatory components (corresponding to constant and electromagnetic- wave magnetic fields, for instance). However, 
the protocol function must satisfy 



t+At 

j f{t')dt' = At, (24) 
t 

and therefore it cannot be purely oscillatory; it must have a constant or other contribution to integrate to a nonzero 
value in accordance with (24). 

The possibility of the modification (23) follows from the fact that the general relation 



U 



time-ordered 



(25) 



does not actually require time ordering as long as the Hamiltonian commutes with itself at different times. This 
condition is satisfied by (23). Furthermore, if the Hamiltonian can be written as a sum of commuting terms then each 
term can he multiplied by its own protocol function. Interestingly, the Hamiltonian of the paramagnetic-resonance 
NOT gatellJ is not of this form. It contains a constant part and an oscillatory part but they do not commute. Note 
that the term proportional to £ in (21) commutes with the rest of that Hamiltonian. The terms proportional to cos 7 
and sin7 do not commute with each other. Rather, they anticommute, in (21), as such terms do in (12). 



6. THE THREE-SPIN XOR GATE 

Thus far we learned that extending the number of spins (qubits, two-state systems) involved in the NOT gate from 
one to two produced an interaction Hamiltonian family (21) with structure that is quite new and unfortunately not 
symmetric in terms of what we are used to in solid-state niagiietic interactions. We will now consider a three-spin 
system: a quantum- XOR gate (which can also be realizedllStjO with two spins). This choice is dictated by the 
fact that we can obtain analytical results and address a new issue that was not there for one- or two-spin systems: 
whether this quantum gate function can be accomplished with two-spin interactions. 

We note that if a quantum logic operation is allowed to be decomposed into a sequence.|Of pUnlimited number of 
universal one- and two-spin gates then one can always reduce it to two-spin interactions.Erl3llj'a Here, however, we 
are interested in one-shot gates for which the external control involves the overall system Hamiltonian, over a single 
time interval Ai. The possibility of using solely two-spin interactions will actually depend on the logical function 
and for more complicated systems it has to be explored by numerical studies. We note also that the issue of having 
the interactions short-range (e.g., nearest-neighbor) does not really arise for few-spin systems although it will be an 
important design criterion as the number of spins (qubits) involved increases. Short-range two-particle interactions 
are much better studied and accessible to experimental probe than multi-particle interactions. 

We denote by A, B, C the three two-state systems, i.e., three spins (qubits). The transformation must be specified 
for those initial states of the input spins A and B, at time that are one of the basis states \AB) = |11), |10), |01), 
or |00), where 1 and denote the eigenstates of the z-components of the spin operators. Here 1 refers to the up state 
and refers to the down state; we use this notation for consistency with the classical bit notion. The initial state of 
G is not specified. We would like to have a quantum evolution that mimics the XOR function: 



A B output 
1 1 

1 1 (26) 
1 1 


Here the output is at time t + At. One way to accomplish this is to produce the output in A or B, i.e., work with 
a two-spin system where the ipflUi and output are the same. The Hamiltonian for such a system is not unique. 
Explicit examples can be foundllJ^LHl where XOR was obtained as a sub-result of the controUed-NOT gate operation. 
In the case of two spins involved, the interactions can be single- and two-spin only. 

Here we require that the XOR result be put in C at time t + At. The final states of A and B, as well as the 
phase of C are arbitrary. In fact, there are many different unitary transformations, C/, that correspond to the desired 
evolution in the eight-state space with the basis \ABC) = |111), |110), |101), |100), |011), |010), |001), |000). The 
choice of the transformation determines what happens when the initial state is a superposition of the reference states, 
what are the phases in the output, etc. 



Let us consider first the following Hamiltoniai 




H = yV2azA<JyB + ViazBCTyC ~ CryB(JxCj ■ (27) 

It is written here in terms of the spin components; the subscripts A, B, C denote the spins. In the eight-state basis 
specified earlier, its matrix can be obtained by direct product of the Pauli matrices and unit 2x2 matrices 2. For 
instance, the first interaction term is proportional to azA (X" cTyB ® , This Hamiltonian involves only two-spin- 
component interactions. In fact, in this particular example A and C only interact with B. 



One can show that the Hamihonian (27) corresponds to the XOR result in C a.t t + At provided A and B were 
in one of the aUowed superpositions of the appropriate-binary states at t (we refer to superposition here because 
C is arbitrary at t). There are two ways to verify this.EllEj Firstly, one can diagonalize H and then calculate the 
evolution matrix U in the diagonal representation by using the general relation (7), valid for Hamiltonians which are 
constant during the time interval At, and then reverse the diagonalizing transformation. 

The second, more general approach presented here is to design a whole family of two-spin- interaction Hamiltonians 
of which the form (27) is but a special case, by analyzing generally a family of 8 x 8 unitary matrices corresponding 
to the three-spin XOR gate. This program is carried out in Sections 7, 8. 



7. THE STRUCTURE OF THE XOR UNITARY MATRIX AND HAMILTONIAN 

We require any linear combination of the states IHl) and |110) to evolve into a linear combination of |110), |100), 
1 010), and |000); compare the underlined quantum numbers with the first entry in (26), with similar rules for the 
other three entries in (26). In the matrix notation, and in the standard basis \ABC) — |111), |110), |101), |100), 
|011), |010), |001), |000), the most general XOR evolution operator corresponding to the Boolean function (26), with 
the output in C, is, therefore. 
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(28) 



The condition of unitarity, UU^ — 1, reduces the number of independent parameters but they are still too numerous 
for the problem to be manageable analytically; we are going to consider a subset of operators of this form. 

From our earlier discussion we know that one way to reduce the number of parameters and ensure unitarity is 
to keep a single phase factor in each column and row of the matrix. Some amount of lucky guessing is involved 
in finding an analytically tractable parametrization. Thus, we choose a form which is diagonal in the states of the 
A- spin, 

^4x4 04x4 
04x4 ^^4x4 



u = 



(29) 



Note that the input spins A and B are not treated symmetrically. Here O4X4 denotes the 4x4 matrix of zeros. The 
4x4 matrices V and W are parametrized as follows: 
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(30) 



(31) 



This choice of an 8-parameter unitary matrix U, see (29), was made because it has the structure 



2U={1 + a,A) V+il- (t,a) W ^V + W + a,A{V - W) , (32) 

where V and W are operators in the space of B and C . Since U was chosen diagonal in the space of A, the 
Hamiltonian, H , will have a similar structure, 



2H = P + Q + a,A{P-Q), 



(33) 



with the appropriate Hamiltonians P and Q in the {B (g) C)-space. In order to avoid three-spin interactions, P — Q 
must be Unear in the Pauli matrices. We also try to avoid single-spin (external-field) interactions. Thus, P + Q must 
contain only terms of the second order in the spin components while P — Q must contain only terms of the first order 
in the spin components. This suggests avoiding putting phase factors on the diagonal, which would lead to matrices 
similar to those encountered in extended-NOT-gate calculations that are known to be of a structure undesirable here: 
they contain a mixture of first-order and second-order terms. The off-diagonal choices remaining are considerably 
limited: the forms (30)-(31) are nearly unique. 

In summary, while the arguments arc admittedly vague and they do involve a certain level of guess, trial and 
error, the presented parametrization offers a good chance that with further restrictions on the parameters a two-spin 
interaction Hamiltonian can be obtained. As will be seen shortly, five conditions are imposed so that the resulting 
Hamiltonian depends on three (real) parameters. 

Let us define 



a + p + j + S p + LJ + ^ + r] 
M = z , = z , (34) 



and introduce the reduced operators p and q, 



Then (7) yields 



V = exp(ip) and W = exp(ig) . 



(36) 



Next, we diagonalize V and W: we calculate their eigenvalues and also the matrices of their normalized eigen- 
vectors, in order to transform to the diagonal representations. 

Specifically, the eigenvalues of V are e*'', ie*^, — e*^, — ie*^. The appropriate eigenvalues of p then follow form 
(36) as fj,, IJ,+ ^TT, H + TT, n+ Itt. Arbitrary multiples of 27r can be added to these choices. However, there are certain 
nonrigorous arguments for generally keeping the spread of eigenvalues of the Hamiltonian as small as possible. Thus, 
we choose the simplest expressions. The eigenvalues of q are determined identically, with jj replaced by u throughout. 

The next step is to apply the inverse of the diagonalizing transformations for V and W to the diagonal 4x4 
matrices for, respectively, p and q. The latter contain the eigenvalues of p and q as the diagonal elements. The 
results are the matrix forms of the operators p and q in the original representation: 



/ |m + 3 

-(1 -z)e*("-^) 
-(l + i)e^(''-'*) 

Y _g»(a;+7-2/:j) 



-(H-i)e'(''-") 

_eH/3+7-2M) 

-(1 - i)e'(T-'') 



-(1 - i)e*(''-^) 

_gi(2A<-/3-7) 

-(l-l-i)e'(''-'') 



-(1 + i)e*(^-'') 
_(l_i)e^(/3-M) 

^M + 3 



(37) 



4 

-q 

TT 



/ 1^ + 3 
-(1-Fi)e*(''-'') 
_(l_i)e*(?~'^) 



v 



ci(2i'-p-uj) _ 



(1 - i)e'^P-''^ 
^ + 3 

_gi(2iy-w-r;) 

(1 -h i)e*(''-'^) 



-(l+i)e 



_^i{uj+ri-2v) 

-(1 - «)e*(''-^) 



_gj(/3+W-2l/) 

-(l + i)e^(''-'') 



(38) 



8. THE TWO-SPIN-INTERACTION XOR HAMILTONIAN 



Thus far we decreased the number of independent parameters in the general unitary transformation and chose it to 
be diagonal in the A-space. We now refine our design of the Hamiltonian to favor interaction of the second order 
in the Pauli matrices. First, we note that both P and Q are constant-diagonal matrices. Therefore, in terms of the 
Pauli matrices both their sum and difference in (33) will involve constant terms. These are undesirable because in 
(Jza{P ^Q) they lead to terms of order one (instead of the desired two), in H, while in P + Q they lead to an additive 
constant in H which only affects the overall phase of the unitary transformation and is of no interest otherwise. 
Therefore, we put 



M = z/ = --7r, (39) 

in order to nullify these diagonal elements in both P and Q. 

Let us now focus our attention on P — Q which, by (39), is now a matrix with zero diagonal. We must impose 
conditions on the parameters to make P — Q oi order exactly one in the Pauli matrices. We note, however, that due 
to zero-diagonal, it cannot contain cr^ terms. The general form liner in Ux-,Oy is 



P-Q=Te 





X* 



X 




+ 





y* 



( ° 


X 


Y 





X* 








Y 


y* 








X 


\ 


y* 


X* 






(40) 



where the stars denote complex conjugation, X and Y are arbitrary (complex) numbers, and X stands for the unit 
matrix as before. Thus we require that P — Q be of the form suggested by (40). This imposes several conditions: 
two above-diagonal elements of the difference must be equal to zero while the remaining four elements must be equal 
pairwise. One can show that these conditions are satisfied provided a, (3, 7 are kept as three independent (real) 
parameters while the remaining angles are given by 



5 = —3tt — a — /3 — 7 , 

p=-n + p, 
LO = —2-K — a — /3 — 7 , 

C = -7r + 7, 
T] = -K + a. 



(41) 
(42) 
(43) 
(44) 
(45) 



Note that (39) is built into (41)-(45). Given this choice, it turns out that P + Q contains only two-spin interactions. 
We have no simple explanation of this property (of the absence of first-order terms in P -|- Q). It is probably related 
to the fact that the structure pattern of the original matrices V and W is quite similar even though the precise 
positioning of nonzero elements is different. We find 



P + Q 



4At 



I 







g«7 _ g«(a-l-/3+7) \/2e?-'^^^'~'^ 
\ \/2e^^"'^'*^ g»("+/3-l-7) _ gi' 



-i(a+[j+^) _ g-i7 





-he 



-i{a+^) 



(46) 



V2Trnil -e*"-he-*^ e-i(a+/3+7) + g-' 

_gi(«-K/3-K7) _ ei7 _g'" + e-i/3 Q 



Finally, we expand these matrices in terms of products of the Pauli matrices and collect terms according to (33) 
to get 



+ [cos{a + 'y) - cos{(3 + 'y)]axB(TyC + [cos(a + 7) + cos(/3 + 7)] (TyBCTxC + [sin(a + 7) - sin(/3 + 7)] crystTyc [ • (48) 



Note that (27) corresponds to the parameter choice a = /3 = 7 = 0. The Hamiltonian (48) describes the three-spin 
XOR. For arbitrary values of a, (3, 7, the interactions involved are all two-spin as desired. The result, however, is not 
symmetric in any obvious way. It seems to correspond to complicated tensor interactions involving expressions of 
order two in the components of the three spins involved. No rotational or other symmetry in the three-component 
spin space, or planar symmetry, or uniaxial coupling, are apparent. All these would correspond to the familiar 
Hcisonbcrg, XY, and Ising couplings in condensed matter physics. 

Thus, in order to realize interaction (48) in materials, a rather anisotropic medium with highly nonsymmetric 
tensorial magnetic interactions will be required. In this respect our analytical attempt to design a multi-spin quantum 
gate cannot be considered really successful. While wc learned several general principles, the imappcaling form of the 
Hamiltonians (21) and (48) suggests that different roots to the derivation of Hamiltonians should be also explored. 
One could start with the more conventional magnetic interactions, isotropic (Heisenberg) , planar (XY), uniaxial 
(Ising), and explore general-parameter Hamiltonians, adjusting the coupling parameters numerically in search of 
those values for which useful Boolean gate operations are carried out. Another approach would be to use, in input 
and output, quantum states more general than the logic products of the qubit-basis-states which are presently 
favored because of the correspondence with classical compiitcrs. Then the design of Hamiltonians may be easier but 
there is a trade-off. These devices will have to be incorporated into larger computers or interact with classical (i.e., 
dissipative) computer components. Radically new programming ideas will be needed to use transformations that 
connect superpositions of qubit states to carry out computations. 
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